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ABSTRACT 

We compare H a emission profiles of 12 dynamically confirmed black holes (BHs) and 2 neutron 
star X-ray transients (SXTs) in quiescence with those of a sample of 43 Cataclysmic Variables 
(CVs), also quiescent. The full-width-half maximum {FWHM) of the H a line in SXTs is tightly 
correlated with the velocity semi-amplitude of the donor star K 2 = 0.233(13) FWHM. This 
new correlation, when combined with orbital periods (i.e. through photometric light curves), 
opens the possibility to estimate compact object mass functions from single integration, low- 
resolution spectroscopy. On the other hand, CVs above the period gap are found to follow 
a flatter correlation, a likely consequence of their larger mass ratios. We also find that the 
FWHM traces the disc velocity at ~42% Rli, independently of binary mass ratio. In addition, 
for a given FWHM, BHs tend to have lower EWs than CVs. This might be explained by the 
fact that CVs must be seen at higher inclinations to mimic the same projected disc velocities as 
BH SXTs. For the same reason CVs with FWHM > 1500 km s -1 are mostly eclipsing while 
none of our sample BHs are. Further, we show that there is a vacant/unoccupied region for 
CVs in the FWHM — EW plane defined by FWHM > 2568 y/(l — ( 9/EW ) 2 ) (krn/s). Both 
the FWHM — K 2 correlation and the FWHM — EW plane can be exploited, together with 
photometric light curves, to efficiently discover quiescent BHs in deep H a surveys of the Galactic 
Plane. 


Subject headings: accretion, accretion disks - binaries: close - Stars: black holes, neutron stars, dwarf 
novae, novae, cataclysmic variables - X-rays: stars 


1. Introduction 


Black holes (BH) provide unique laboratories to 
study key astrophysical phenomena such as accre¬ 
tion, the production of relativistic outflows and 


gamma ray bursts (e.g. Fender & Gallo 2014 


Brown et al. 2000). BHs are usually discov¬ 


ered in transient X-ray binaries (SXTs) thanks to 
their dramatic, large-amplitude X-ray outburst^] 
These are caused by thermal-viscous instabilities 
in accretion discs which are (slowly) fed by matter 


transferred from a (donor) companion star (Lasota 


2001). The distribution of BH masses and binary 


periods provide fundamental constrains to models 
of supernovae explosions and compact binary evo- 


1 By contrast, persistent X-ray binaries tend to harbour neu¬ 
tron stars. 


lution. Unfortunately, our understanding of the 
physics involved in these processes is still patchy 
because observational data are severely limited by 
small samples. 


The current sample of ~50 BH SXTs, with only 
17 dynamical confirmations, is the tip of the ice¬ 
berg of a hidden population of hibernating BHs 
(see Casares & Jonker |2014 1. The size of the over¬ 
all population is highly uncertain. Extrapolation 
of the number of SXTs discovered in the X-ray era 
suggests that several thousand dormant BHs are 


awaiting to be discovered (van den Heuvel 1992 


Tanaka fe Shibazaki||1996~ Romani||1998 ). On the 

other hand, modern population-synthesis models 
predict an even larger population of ~ 10 4 — 10 5 


SXTs (Pfahl et al. 2003 Yungelson et al. 2006). 


Observational constraints are likely biased low be¬ 
cause they suffer from incompleteness and neglect 
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a hitherto unexpected population of long period 
SXTs with very long duty cycles or even supressed 
outburst activity (cf. Ritter & King 2002; Menou 


et al. 19991. Furthermore, there is mounting ev¬ 


idence for the existence of a population of X-ray 


obscured or intrinsically faint BH SXTs (Corral- 

Santana et al. 

2013 Armas Padilla et al. 

2014). 


It has been shown that the latter could be mem¬ 
bers of a sizeable population of short period BH 


transients with low outburst luminosities (Mac- 


carone & Patruno 2013). And finally, the recent 


discovery of the first Be/BH binary indicates that 
some BHs might be accreting ’’silently” from the 
slowly outflow winds of rapidly rotating Be-type 
stars (Casares et al. 20141. While binary popu¬ 
lation models predict a very modest number of 


Be/BH binaries in the Milky Way (Belczynski & 

Ziolkowskij 

2009 

Grudzinska et al. 

2015 

) empiri- 


cal constraints, with only one detection yet, are 
very loose. Altogether, the only way to make sig¬ 
nificant progress in our knowledge of the Galactic 
population of BHs requires the discovery of a large 
sample of quiescent SXTs and this demands a new 
research methodology. 

Quiescent BH SXTs are particularly difficult to 
find since they are relatively faint across the elec¬ 
tromagnetic spectrum in this state. By their very 
nature, quiescent states are characterised by ex¬ 
tremely low accretion levels (< 10~ 9 M 0 yr _1 ). 
The inner disc is truncated in an advected flow 
and, further, the lack of a solid surface (an exclu¬ 
sive signature of BH) results in extremely weak 


X-ray, UV and radio luminosities ( 

Narayan & Me- 

Clintock 

2008; Miller-Jones et al. 2011 

). On the 


other hand, the optical spectrum is dominated by 
the low-mass donor star, with superposed emis¬ 
sion lines from the accretion disc gas. The lack 
of a hard radiation field implies that only emis¬ 
sion from neutral H and He is detected. Crucially, 
the strongest emission line is H a , but several other 
Galactic populations are also strong H a emitters, 
such as Cataclysmic Variables (CVs), Symbiotic 
binaries, flare stars, reddened Be stars, T Tauri 
and other classes of young stellar objects. As a 
consequence, H Q surveys of the Galactic plane are 
vastly outnumbered by these populations of con¬ 
taminating H q emitters. Attempts to clear out the 


sample using color selection cuts (Corral-Santana 


et al. 2008) or X-ray diagnostics (Jonker et al. 


2011) still have to prove their effectiveness. Inter¬ 


estingly, new generation radio surveys (SKA and 
its pathfinders) may offer an alternative route to 
detect quiescent BHs given the increase in radio- 


to-X-ray flux ratios at very low luminosities (Mac- 

12005 1. 


carone 


Here we present the discovery of a correlation 
between the width of the H Q line and the pro¬ 
jected velocity of the companion star in quiescent 
SXTs. This property can be applied, in combina¬ 
tion with photometric orbital periods, to gather 
compact object mass functions and flag new po¬ 
tential BHs. This strategy can be turned into a 
novel technique to unveil hibernating BHs, tech¬ 
nique which appears much more efficient than tra¬ 
ditional methods based on time-consuming spec¬ 
troscopic monitoring of new X-ray novae. 


2. The sample 
2.1. X-ray Transients 

We have assembled a spectroscopic database 
of dynamically confirmed BH SXTs with H a 
emission (see Table 1). Most spectra were col¬ 
lected by us (V404 Cyg, BW Cir, N Mus 91, 
GS 2000+25, A0620-00, XTE J1650-500, XTE 
J1859+226, GRO J0422+320, XTE ,11118+480) 
and have been presented in several publications, 
while others were kindly provided by J. Orosz 
(XTE J1550-564), R. Remillard (N Oph 77) and 
A. Fillipenko (N. Vel 93). In addition, three new 
unpublished spectra of GRO J0422+320 were ob¬ 
tained on the night of 28 Jan 2009 with ALFOSC 
on the 2.5m Nordic Optical Telescope (NOT) at 
the Observatorio del Roque de los Muchachos. In¬ 
tegration time was set to 2400s and we used grisrn 
#4 and 1” slit to cover the range AA3820-9140 at 
360 km s -1 resolution. Seventy unpublished H a 
spectra of V404 Cyg were also collected between 
1994-2005 with the 2.5m Isaac Newton Telescope 
(INT) and the 4.2m William Herschel Telescope 
(WHT) at resolutions 36—134 km s -1 . Five ad¬ 
ditional spectra of V404 Cyg were obtained on 
the nights of 5 and 7-8 July 2008 at 200 km s -1 
resolution with the 2.1m telescope at the Observa¬ 
torio San Pedro Martir (SPM). Finally the sample 
also includes spectra of the two only neutron star 
SXTs with a reported radial velocity curve of the 
companion star, Cen X-4 and XTE J2123-058. 
Table 1 provides the main observational details 
and associated references for every system. 
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2.2. Cataclysmic Variables 


3. The FWHM — K 2 correlation 


We also collected a database of H a spectra of 
Cataclysmic Variables (CVs) in quiescence includ¬ 
ing 41 dwarf novae and two intermediate polars, 
also classified as classical novae (GK Per and DQ 
Her). Most of the spectra were acquired during 
several observing campaigns and Service nights 
performed with the WHT in 1992, 1993, 1998 and 
2006, the INT in 1992, 2001, 2008 and 2009, the 
NOT telescope in 2008, 2009 and 2012, the 2.2m 
telescope at Calar Alto Observatory in 1995 and 
the 2.1m telescope at SPM in 2008. The WHT 
spectra were obtained with ISIS and the R1200R 
grating which delivers 40 km s _1 resolution at 
H a . The INT spectra were obtained with the 
IDS spectrograph and gratings R150V, R300V, 
R900V and R1200R covering the H Q region at 50- 
320 km s -1 resolution. For the NOT campaigns 
we used ALFOSC and grism #7 with different 
slit widths providing 235—310 km s -1 resolu¬ 
tion. The Calar Alto data were collected with 
the CASSPEC spectrograph, the f3 camera and 
grating #11 in second order which delivered 27 km 
s -1 resolution. The SPM spectra were obtained 
with the Boiler & Chivens spectrograph and a 
600 1/mm grating to yield 192 km s -1 resolution. 
Spectra of 13 CVs (BV Cen, EY Cyg, DX And, 
HS0218+3229, AH Her, HS2325+8205, SDSS 
J 100658.40+233724, U Gem, CTCV J1300-3052, 
OY Car, V2051 Oph, SDSS 103533.02+055158 
and SDSS J143317.78+101123.3) were kindly pro¬ 
vided by different authors while H a parameters for 
other 15 CVs were compiled from literature. In 
addition, spectra of WZ Sge were retrived from the 
ING archive and reduced by us. The final collec¬ 
tion of CV spectra were employed as a test sample 
for comparison with the H Q properties of SXTs. 
All the CVs have secured velocity amplitudes (iv 2 ) 
of their companion stars either through radial ve¬ 
locity studies or eclipse light curve solutions (for 
short period binaries where the companion star is 
undetected). Table 2 summarizes the main obser¬ 
vational details of the CV database. The database 
is intended to be a representative sample of qui¬ 
escent CVs with available K 2 determinations al¬ 
though we warn about possible selection effects. 
These will be addressed in the following Section. 


Full-width-half-maximum {FWHM) values 
were obtained from gaussian fits to individual H a 
profiles in every SXT. The fitted model consists of 
a constant plus a Gaussian function. Continuum 
rectified spectra were fitted in a window of ±10000 
km s^ 1 , centered on the H a line after masking the 
neighboring Hel line at A6678. We adopted 1 -a 
formal errors on the fitted parameter as derived 
through x 2 minimization. Fig. 1 displays some 
fit examples to average line profiles covering the 
entire range of line widths displayed by our data. 
It is clear from the figure that a simple Gaus¬ 
sian does not provide an accurate description of 
very broad profiles with large double peak separa¬ 
tions. However, we find that the FWHM given 
by the Gaussian model is within 10% of other 
width parameters obtained from more sophisti¬ 
cated double-Gaussian models. And, more impor¬ 
tantly, it is far more robust since double-Gaussian 
models can easily fail when fitting low signal-to- 
noise profiles. In addition to the FWHM , we also 
extracted equivalent widths (EW) by integrating 
the H„ flux in individual spectra, after continuum 
normalization. 


Our 266 spectra of V404 Cyg span over 20 
yr and, therefore, present the most complete 
database yet for the analysis of the secular evo¬ 
lution of accretion discs in quiescent BHs. Fig. 
2 presents the evolution of both line parameters, 
FWHM and EW, in V404 Cyg. The data points 
have been folded into 50 day bins, except for the 
first 30 days where 1 day bins were used to better 
trace the rapid evolution through the outburst. 
The plot shows a steep rise in FWHM (drop in 
EW) followed by a plateau phase starting ~1300 
days after the peak of the outburst. By taking 
the FWHM as a proxy of the disc radius, Fig. 
2 indicates that the accretion disc shrinks during 
outburst, recovering an equilibrium radius ~3.5 
yr after the maximum. Evidence of accretion disc 
shrinkage following outburst has been reported 
for CVs using eclipse timing and eclipse mapping 
techniques (Srnak 1984 Baptista & Catalan 2001). 
Because we are interested here in comparing av¬ 
erage properties of stable quiescent discs in SXTs 
we decided to trim all data obtained within ~1300 
days after the peak of the outburst. This leaves 
127 useful spectra of V404 Cyg (all since 1993), 94 
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of BW Cir and 3 of GRO J0422+320 (from 2009). 

Table 3 lists the mean FWHM and EW val¬ 
ues for SXTs, where the quoted uncertainties rep¬ 
resent 1 standard deviation in the distribution of 
individual measurements. Therefore, our errors 
mostly reflect time variability in line width and 
normalized flux. Line variability is mainly caused 
by aperiodic flares (Hynes et al.|2002) or long-lived 
disc asymmetries (e.g. hot spots), modulated with 
the orbital period. In the cases of N. Oph 77 and 
Vel 93 only one phase averaged spectrum is avail¬ 
able and, thus, we adopt gfwhm = 0.1 FWHM 
to account for the average 10% variability dis¬ 
played by the remaining binaries. By the same to¬ 
ken, we adopt (Jew = 0.22 EW for objects where 
only one spectrum is available, based on the mean 
standard deviation in the EW s of the remain¬ 
ing systems. Instrumental resolution was sub¬ 
tracted quadratically from every FWHM value 
and, therefore, the H Q widths quoted in Table 3 
are intrinsic. Table 3 also provides fundamen¬ 
tal binary parameters, chiefly the orbital period 
(P), the radial velocity semi-amplitude of the com¬ 
panion star (K 2 ) and the mass of the compact 
star (Mi) and inclination (i), when available, with 
their associated references. 


We also collected FWHMs and EWs from 
Hq, lines in our sample of CV spectra listed in 
Table 2. In the case of 14 eclipsing binaries 
(EM Cyg, EX Dra, HS 2325+8205, DQ Her, 
SDSS J100658.40+233724.4, U Gem, IP Peg, 
CTCV J1300-3052, HT Cas, OY Car, V2051 
Oph, SDSS 103533.02+055158.3, WZ Sge and 
SDSS J143317.78+101123.3) only spectra ob¬ 
tained ±0.05 phases away from the eclipse min¬ 
imum were considered. As already mentioned, 
FWHMs and EWs values were obtained from 
the literature for 15 CVs. For these cases we 
adopt a 7% error in FWHM and 14% error in 
EW, derived from the mean variability measured 
in the other 28 CVs. These values are listed in Ta¬ 
ble 4, together with determinations of the orbital 
period and K 2 velocities. 

In Figure 3 we compare K 2 versus FWHM and 
it can be seen that these quantities are tightly cor¬ 
related in SXTs, with a Pearson correlation coef¬ 
ficient r = 0.99. A linear fit yields the following 
relation 


where both quantities are given in km s _1 . A sec¬ 
ond order polynomial fit is not justified since the 
constant coefficient is consistent with zero at ler. 
In order to assess the error in K 2 estimated from 
the correlation we computed the difference with 
the true (dynamical) K 2 for our 14 SXTs. The 
values are found to follow a Gaussian distribution 
with standard deviation of 22 km s -1 . We there¬ 
fore conclude that robust estimates of the K 2 ve¬ 
locity can be obtained from the width of the H a 
line in SXTs. The uncertainty in the coefficient of 
the correlation was estimated through a Monte- 
Carlo simulation of 10 4 events, assuming that the 
difference between the model and true I\ 2 values 
follow a Gaussian distribution with a = 22 km 
s- 1 . 

The FWHM — K 2 correlation is expected from 
basic equations. Assuming that FWHM is de¬ 
termined by gas velocity at a characteristic disc 
radius f?w 


/ FWHM \ 2 GMi . 2 

HH =iiw ra ‘ <2) 

where Mi is the mass of the accreting star and 
i the binary inclination. On the other hand, the 
companion’s velocity is given by 


Rk 2 = 


GMl 


a {Mi + M 2 ) 


■ 2 • 
sm 1 


(3) 


with M 2 the mass of the companion star. There¬ 
fore, 


( I < 2 

l FWHM 


Rw 

4 a (1 + g) 


(4) 


where q = M 2 /Mi is the mass ratio and a the 
binary separation. If we now assume Ryj = oPli 


(with a < 1) and use Eggleton’s relation (Eggleton 
1983) to remove Rli/cl then 


l<2 

FWHM 


2 


where 


(5) 


m = 


0.49 (l + g)" 


0.6 + q 2 / 3 In (l + g -1 / 3 ) 


( 6 ) 


K 2 = 0.233(13) FWHM 


Cl) 


4 


In the domain of SXTs, with q ~ 0.05 — 0.15, 
the dependence of Eq. 5 with q is modest because 














\Jf(q) varies between 0.77—0.69. For typical BH 
SXTs, with q = 0.1, K 2 /FWHM = 0.36^0 and 
hence, the empirical correlation, described by eq. 
1, implies that the FWHM of the H„ line in BH 
SXTs traces the disc velocity at about 41% Rr/\ ■ 
And given that ~ 80% of the line flux is contained 
within one FWHM and that typical quiescent 
disc radii reach 0.5 — 0.6i?£i ( Marsh et al.|1994 


Casares et al.||1995b) we observe that the bulk of 


the H a emission arises from the outermost regions 
of the accretion disc. 

Equation 1 stems from the fact that the Ke- 
plerian velocity field in the H a disc provides a 
fundamental dimension scale of the binary. Sim¬ 
ilarly to the rotational broadening of the donor 
star V sin i (see Wade &; Horne 1988), the mean 
accretion disc velocity (traced by the FWHM) 
scales with the donor star’s velocity. Since both 
are projected velocities along the line of sight, the 
dependence on inclination cancels out. However, 
unlike V sini, this new correlation is only weakly 
dependent on binary mass ratio and hence results 
in a very tight linear regression. We also note 
that eq. 1 is more robust than a former empirical 
relation between K 2 and the double peak sepa¬ 
ration (Orosz et al. 1994 Orosz & Bailyn 19951 
because the latter traces the tidal (outer disc) ra¬ 
dius which can be strongly affected by disc asym¬ 
metries, the presence of S-wave distortions driven 
by hot-spots and anomalously low (sub-Keplerian) 
velocities (e.g. North et al. (2002)). 

At this point it is instructive to compare how 
the CVs distribute in the FWHM — K 2 plane. 
But prior to this it is important to asses the im¬ 
pact of possible selection effects. To start with, 
since the companion star is mostly detected above 
the period gap, our reference sample is, in princi¬ 
ple, biased toward large mass ratios. The effect 
is further exacerbated by the fact that the CV 
population in the 3-4 h period range is dominated 
by SW Sex and Novalike stars, systems in perma¬ 
nent outburst (Rodrfguez-Gil et al. [2007). This 
explains the paucity of CVs with periods < 0.17 
d in our list. To compensate for the deficit of bi¬ 
naries with small q values we made an effort to 
incorporate short period CVs. These are, how¬ 
ever, strongly skewed towards high inclinations be¬ 
cause only eclipsing CVs can yield reliable I\ 2 val¬ 
ues (through light curve modeling) when the com¬ 
panion star is not detected. Fig. 4 presents his¬ 


tograms of the distribution of orbital periods, mass 
ratios and inclinations in our CV sample, clearly 
depicting these selection effects. Mass ratios and 
inclination values were compiled from the refer¬ 
ences listed in Table 4 and Ritter & Kolb (20031. 
In the figure we make a distinction between CVs 
above and below/within the period gap. 

Figure 4 shows that the distribution of mass ra¬ 
tios is clearly bimodal and can be described by two 
Gaussians: 32 long-period CVs cluster at q = 0.63 
with a = 0.2 while 11 short-period CVs define a 
much narrow peak centered at q = 0.12 and with 
a = 0.07. We note that this is overall similar to 
the distribution of mass ratios obtained from the 
entire CV data available in the Ritter’s catalogue 
(Ritter & Kolb 2003). The figure also highlights 
the fact that our short period CVs are strongly bi¬ 
ased toward high inclinations. By constrast, there 
seems to be no significant bias in the inclination of 
CVs above the gap. In view of this, we find justi¬ 
fied to distinguish hereafter between short-period 
CVs (i.e. below/within the gap) and long-period 
CVs (above the gap), with the latter not being 
strongly affected by selection effects. 


The bottom panel in Fig. 3 presents the loca¬ 
tion of the CVs in the FWHM — K 2 plane. The 
figure shows that long-period CVs display a sim¬ 
ilar regression to that found for SXTs albeit flat¬ 
ter, i.e. for a given /\ 2 , H a lines are systemat¬ 
ically broader than in SXTs. A linear fit yields 
K 2 = 0.169(16) FWHM , relation which could be 
used to infer K 2 velocities for quiescent CVs above 
the period gap. We attribute the flatter slope of 
the CV correlation to their comparatively larger 
mass ratios, which leads to smaller R^i and thus 
smaller disc radii (in binary separation units). 

If we now bring q=0.63 (i.e. the peak in the q 
distribution of long-period CVs) into eqs. 5-6 and 
set K 2 /FWHM = 0.169 from the empirical fit 
one obtains a = 0.43, in excellent agreement with 
what was found for SXTs. This implies that the 
FWHM of the H a line is always formed at about 
42% Rli, irrespectively of the binary mass ratio. 
Despite the large spread in mass ratios the CV 
correlation appears quite narrow, a consequence of 
the very weak dependence of K2/FWHM with q 
for large q values (see dotted lines in the bottom 
panel of Fig. 3). 

The group of short-period CVs, on the other 
hand, concentrate at large FWHM values be- 
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cause of their high inclinations. They are seen to 
depart from the trend defined by the long-period 
CVs, approaching the SXT correlation, a result of 
their small q-values. Although the short-period 
CVs represent a very biased sample, they are of 
particular interest because they define the upper 
limit in the FWHM distribution of the CV pop¬ 
ulation. 

4. The EW - FWHM plane 

In a given system, the FWHM of any emis¬ 
sion line depends on the binary inclination, orbital 
period and the mass of the compact object. By 
bringing K 2 = (27rasin*)/P(l + q ) into eq. 5 and 
using Kepler’s Third Law we find 


FWHM oc 


sin* 

yja (1 + q) 4/3 f(q) 



(7) 


Here, the dependence on mass ratio is extremely 
weak, with i/(l + g) 4 / 3 /(g) varying in the range 
0.80 to 0.69 for q = 0.05 — 1. Therefore, we can 
safely assume 


FWHM ~ A 


Mi 


1/3 


Sill* 


( 8 ) 


where A is a constant that can be calibrated using 
P and FWHM values listed in Table 3, together 
with dynamical masses and inclinations available 
in the literature for 9 SXTs (also listed in Table 
3). These yield a mean value A = 876 ± 48 km 
s -1 , when M\ and P are expressed in units of Mg 
and days, respectively. 

On the other hand, the EW of the H a line 
depends on the binary geometry. For instance, 


Warner (19861 showed that the EW of H^ in CVs 
increases with inclination because of the reduc¬ 
tion of continuum brightness as the disc is seen 
at large inclinations. Therefore, to a first-order 
approximation one can assume 


EW 


B 

cos i 


(9) 


where the constant B = 9 ± 8 A has been cali¬ 
brated using the distribution of EW s and inclina¬ 
tions listed in Table 3 El 


Fig. 5 displays our sample of SXTs and CVs in 
the EW—FWHM plane. Here we have used open 
triangles to mark eclipsing CVs. BHs are a factor 
~ 10 more massive than white dwarfs or NS and 
thus should possess, on average, wider H„ lines by 
a factor ~ 2.2. Instead, we observe that only one 
BH (XTE J1118+480) stands out clearly in the 
right side of the diagram with FWHM > 2500 
km s -1 . The remaining BHs are mixed up with 
CVs at lower FWHM values because they ei¬ 
ther have long orbital periods (i.e. cases of V404 
Cyg and BW Cir with FWHM ~ 1000 km s" 1 ) 
or are viewed at lower inclinations. Fortunately, 
they cluster in the central part of the diagram 


gion populated by eclipsing CVs. Since the lat¬ 
ter are easily detected through deep (~2-3 mag) 
eclipses we conclude that non-eclipsing binaries 
with FWHM > 1500 km s -1 are good candi¬ 
dates to host BHs. Incidentally, a population of 
(short-period) eclipsing BHs with very wide H a 
lines would be expected in the right part of the 
diagram. Scaling from the eclipsing CVs we pre¬ 
dict them to show a factor ~ 2.2 larger widths 
i.e. FWHM ~ 4200 km s -1 . As a matter of 
fact, the transient X-ray binary Swift J1357-0933 
has been proposed as an extreme inclination BH, 
although no dynamical solution is yet available 
(Corral-Santana et al. 2013). We have recently 
obtained 4 H a spectra of Swift J1357-0933 with 
OSIRIS and the R300R grism on the 10.4m Gran 
Telescopio Canarias (GTC) on 29-30 June 2013, 
from which we measure FWHM = 4085 ± 328 
km s -1 and EW = 131.9 ± 14.5. These values 
are indeed consistent with a short period BH seen 
at very high inclination. The position of Swift 
J1357.2-0933 is indicated by an asterisk in the di¬ 
agram. 


Using eqs. 8-9 one can define regions of con¬ 
stant Mi/P and inclination in the EW — FWHM 
plane. These are marked in the bottom panel of 
Fig. 5 using dashed and dotted lines respectively. 
We stress, however, that the quoted inclinations 
must be considered as mere indicative given the 
crude approximation involved in eq. 9. The fig¬ 
ure suggests (with all the caveats associated to low 
number statistics) that, for a given FWHM , BHs 


in literature for a subset of 25 CVs from our Table 4 we find 
A = 897 zb 46 km s —1 and B = 7 =b 5 A, in good agreement 
with the calibration obtained using SXTs. 


2 If we use instead reliable masses and inclinations reported 
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tend to have lower EWs than CVs. This could be 
explained because, due to their shallower potential 
wells, CVs must be seen at higher inclinations to 
mimic the same projected disc velocities as BHs 
and, therefore, their H a fluxes are less diluted by 
the accretion disc continuum. It is interesting to 
note the position of GRO J0422+320 in the up¬ 
per left side of the diagram. While the mass of its 
BH is quite uncertain (see Casares & Jon ke~r| 2014) 
both the large EW and low M-\ /P factor suggest 
it hosts a low-mass BH, in line with the results 
of Casares et al. (1995a) and |Gelino & Harrison 
([2003F! 


Finally, we can tentatively define a forbidden 
region for CVs in the FWHM — EW plane by 
taking extreme parameters i.e. M\ >1.4 M 0 and 
P < 80 min, the period minimum spike observed 


in the distribution of CV periods (Gansicke et al. 
2009). This yields Mi/P > 25.2 M 0 /d and thus 
FWHM > 2568^(1 - (9 /EW) 2 ), limit which is 
marked by a green solid line in the plot. As a test 
we have examined a random sample of 236 dwarf 
novae selected from Sloan DR7. Sloan CVs show 
the same trend as seen in Fig. 5 i.e. they populate 
the region leftwards of the green line, with a large 
spread towards high EWs up to ^350 A. 


5. Discussion: new strategies to detect 
dormant BHs 

To make progress in our understanding of the 
formation and evolution of Galactic BHs it is 
essential to discover a large sample of secured 
(dynamically confirmed) BHs. This new sample 
would allow us to constrain the number density, 
the orbital period distribution and, ultimately, the 
BH mass spectrum. Only then we will be able to 
address fundamental questions such as the role of 
supernova models in shaping the distribution of 


nn masses, a cm 

(Ozel et al. 2010 

rent not tc 

Farr et al. 

>pic m tne community 
2011 Kreidberg et al. 

2012 Belczynski et al.|2012 

). Deep H a surveys of 


the Galactic plane, combined with spectroscopic 
surveys, can efficiently select samples of H Q emit¬ 
ting objects with FWHM > 1000 km s -1 . This 
width cut would allow instant removal of Galactic 
populations of narrow H„ emitters such as plan¬ 
etary nebulae, Be, chromospheric stars, T Tauri 
and other YSOs. Only high inclination CVs, due 
to the large gravitational fields of their accreting 


white dwarfs, are able to produce wide H Q lines. 

Fortunately, as we have shown, the width of 
the H a line in quiescent BH and NS SXTs is 
tightly correlated with the projected velocity of 
the donor star. The FWHM — K 2 relation can, 
therefore, be exploited, together with supplemen¬ 
tary information on orbital periods (e.g. from 
light curve variability), to gather ’’preliminary” 
mass functions ( PMF ) of compact objects from 
single epoch spectroscopy i.e. PMF = 1.3 x 
10 ~ 9 p (FWHM ) 3 , where P is given in days 
and FWHM in km s -1 . K 2 values can be esti¬ 
mated from single integrations rather than expen¬ 
sive time resolved spectroscopy, allowing a search 
for dynamical BHs in much deeper fields and using 
a factor ~4 lower spectral resolution than usually 
employed. Therefore, the novel strategy that we 
propose is clearly much more efficient than stan¬ 
dard spectroscopic techniques, aiming at measur¬ 
ing the orbit of the donor star from the Doppler 
shift of weak absorption lines. And it may be the 
only way to infer mass functions in extremely faint 
BH SXTs i.e. the bulk of the Galactic population. 
We also note that our technique can be easily ex¬ 
ecuted in crowded regions like globular clusters, 
where HST time-series photometry can yield or¬ 
bital periods while ground-based spectroscopy is 
severely limited by seeing conditions. 

Fig. 6 demonstrates how BHs are nicely segre¬ 
gated from CVs and NS using the FWHM — K 2 
correlation. Every BH (except GRO J0422+320) 
is located in the right hand part of the diagram, 
beyond the dotted vertical line. Therefore, tar¬ 
gets with PMF > 2 M 0 are strong candidates 
to host BHs. It should be noted that the plot¬ 
ted CV PMFs are, in most cases, robust upper 
limits to true mass functions because, by adopt¬ 
ing eq. 1 we overestimate their real K 2 values. 
By comparing the PMF values with the true 
(dynamical) mass function of BHs one can esti¬ 
mate the typical uncertainty introduced by the 
FWHM — K 2 correlation. Relative errors in mass 
function are found to follow an approximate Gaus¬ 
sian distribution with <7 = 0.08. Therefore, the 
FWHM — K 2 correlation allows us to estimate 
mass functions with typically ~10% uncertainty. 
As an example, we have applied eq. 1 to Swift 
J1357.2-0933, where our quiescent GTC spectra 
yield FWHM = 4085 ± 328 km s -1 . Extrapolat¬ 
ing from eq. 1 we predict K 2 = 952 ± 93 km s _1 
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which, when combined with P = 0.117 ± 0.013 d, 
(Corral-Santana et al. 2013) leads to a record mass 
function PMF = 10.5+3.3 Mq. We note that this 
figure might even be slightly underestimated since 
the GTC spectra were taken only ~ 900 d. after 
the peak of the outburst (see Sect. 3). In any 
case, the effect would be small compared to our 
errorbar, which is dominated by propagating the 
large FWHM uncertainty (driven by line variabil¬ 
ity) into the FWFIM — I\ 2 relation and the mass 
function equation. 


A big step forward in the exploitation of the 
FWFIM — K 2 relation will likely come from the 
use of imaging techniques. Accurate H Q widths 
can be measured directly through a combination 
of customized narrow-band filters, removing the 
need for any spectroscopy at all. This brings in 
a new observational signature that we coin here 
Photometric Mass Function (PMF, taking advan¬ 
tage from the same acronym before). PMF s open 
a novel concept, i.e. that of weighting mass func¬ 
tions photometrically, and lay the ground for the 
efficient discovery of new hibernating BHs in large 
survey volumes. We are currently working on this 
strategy. 


Finally, it should be borne in mind that there 
is a strong selection effect against detecting high 
inclination BHs in X-ray selected samples. This 
is because flared accretion discs obscure X-rays 
and indeed none of the currently known dynam¬ 
ical BHs has an inclination > 75° (Narayan & 


McClintock 20051. Since a PMF-b&sed sample 
will be selected by H Q widths and not X-ray emis¬ 
sion we expect this strategy will uncover a signif¬ 
icant number of high inclination and even eclips¬ 
ing BHs. These hold the prospect to render the 
most accurate BH masses yet because the relative 
uncertainty in the inclination dominates the mass 
error budget. Clearly, the newly discovered BHs 
will have an strategic impact in the construction 
of the BH mass spectrum. 
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Fig. 1. — Example of Gaussian fits to Ha profiles 
in SXTs. A selection of average spectra, repre¬ 
senting the entire range of FWHM s, is depicted. 
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Time since outburst (days) 


Fig. 2.— Long-term (20 year) evolution of the 
FWHM and EW of the H Q line in V404 Cyg. 
Individual points have been co-added into 50 day 
bins except for the first 9 data points where 1-day 
bins were used. Errorbars reflect variability within 
the bin. 
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Ha FWHM (km/s) 



Ha FWHM (km/s) 


Fig. 3.- Top: the FWHM — K 2 correlation for 
SXTs with the best linear fit. Blue solid circles 
indicate BHs while Ns are marked by open circles 
Bottom: the FWHM — K 2 correlation for CVs. 
Filled triangles indicate CVs above the period gap 
while open triangles those below/within the gap. 
Like in the case of SXTs, long period CVs display 
a tight FWHM — K 2 correlation although with a 
flatter slope. Dotted lines mark theoretical corre¬ 
lations from eq. 5-6 for q = 0 — 1 and a = 0.43. 
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Fig. 4.- Distribution of orbital periods (top 
panel), mass ratios (middle) and inclinations (bot¬ 
tom) for the test sample of CVs. The white his¬ 
togram represents CVs with orbital periods above 
the period gap while the red shaded histogram 
those below/within the gap. The period gap is in¬ 
dicated by dashed vertical lines in the top panel. 

All the CVs from the white histogram possess if 2 
determinations based on radial velocity curves of 
the donor stars while those in the red shaded his¬ 
togram on modeling eclipsing light curves. The 
only exceptions to the latter are the ultra-short 
period system El Psc and the period gap CV SDSS 
J1300-3052. It can be seen that short period CVs 16 
have very low mass ratios and are strongly biased 
toward high inclinations (with the exception of the 
2 CVs just mentioned) . 
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Fig. 5.— Top: SXTs and CVs in the FWHM - 
EW plane. Same symbol code is used as in Fig. 
3, except for open triangles which now indicate 
eclipsing CVs. Bottom panel: dashed vertical lines 
correspond to constant Mi/P or b values (expressed 
in units of Mg/d) while horizontal dotted lines 
indicate various inclinations, in an idealized case 
where EWs are fortshortended by a factor cos* 
due to disc continuum visibility. The green solid 
line defines the maximum FWHM predicted for 
CVs. 
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Fig. 6.— Preliminary mass functions obtained 
from K 2 values derived through the FWHM — K 2 
correlation for SXTs. Same symbol code is used 
as in Fig. 5. 


18 





























Table 1 

Database of X-ray transients 


Object # Spectra Year Resolution (km/s) ref 


Black Holes 


V404 Cyg 

266 

1989-2009 

6-180 

1-8 

BW Cir 

96 

1995-2006 

70-110 

9-10 

XTE J1550 -564 

33 

2001, 2008 

55-165 

11-12 

N. Oph 77 

1 

1993 

180 

13 

N, Mus 91 

29 

1993-1995 

74-90 

14 

GS 2000+25 

25 

1995 

196 

15 

A0620-00 

20 

2000 

7 

16 

XTE J1650-500 

15 

2002 

35 

17 

N Vel 93 

1 

1998 

120 

18 

XTE J1859+226 

10 

2010 

255 

19 

GRO J0422+320 

21 

1994-1995, 2009 

230-630 

8, 20 

XTE J1118+480 

120 

2011 

120 

21 

Neutron Stars 

Cen X-4 

90 

1993-2002 

6-74 

22-24 

XTE J2123-058 

20 

2000 

123-183 

25 
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(6) Hynes et al. (20021; (7) 


Casares et al. (1993 1; (5) jCasares fe Charles ( |1994 ); 


Gonzalez Hernandez et al.|(|2011); (8) This paper; (9) 


Casares et al. 

(2004); (10) 

Casares et al. 

(2009a); (11) Orosz et al. 

(2011) 

Orosz et al. (2002); (13) Remillard et al. (1996|; (14) Casares et al. 

(1997) 


al. (2004); (18 ) |Filippenko et al. (1999); (19) Corral-Santana et al. ( 2011|); (20) 


Casares et al. (1995a); (21) Gonzalez Hernandez et al. (20121;(22) 


Torres et al. 


(2002); (23) D’avanzo et al. (2005); (24) Casares et al. (2007); (25) 
( 2002 ) . 


Casares et al. 
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Table 2 

Database of Cataclysmic Variables. 


Object 

# Spectra 

Year 

Resolution (km/s) 

ref 

GK Per 

8 

1992, 1995, 1998 

25-70 

1-2 

SDSS 2044-04 

- 

2001-2004 

190-210 

3 

BV Cen 

63 

2004 

10 

4 

RX 1951.7+3716 

- 

2001-2004 

190-210 

3 

UY Pup 

- 

2002-2003 

190-210 

5 

EYCyg 

103 

1998-2001, 2004 

15-20 

6 

DX And 

4 

1993 

180 

7 

SY Cnc 

28 

1992, 1998-1999, 2003, 2008 

25-70 

1,8 

AT Ara 

- 

2000 

100 

9 

RU Peg 

2 

2008 

60 

2 

GY Hya 

- 

2001-2004 

190-210 

3 

CH UMa 

16 

2008-2009 

60-320 

2 

V392 Hya 

- 

2001-2004 

190-210 

3 

RY Ser 

- 

2003 

190-210 

5 

HS 0218+3229 

56 

2002 

55-75 

10 

EM Cyg 

3 

2008 

60 

2 

Z Cam 

12 

2009 

60 

2 

SDSS 0813+45 

- 

2002 

190-210 

5 

V426 Oph 

4 

2006 

40 

2 

ss Cyg 

23 

2008 

60-250 

2 

LY UMa 

- 

2000 

190-500 

11 

BF Eri 

5 

2009 

235 

2 

TT Crt 

- 

2000-2003 

190-210 

5 

AH Her 

10 

1980-1981 

60 

12 

EX Dra 

31 

2001 

50 

2 

HS 2325+8205 

42 

2005, 2007 

215 

13 

DQ Her 

4 

2009 

60 

2 

SDSS J100658.40+233724 

68 

2008 

140 

14 

SS Aur 

28 

2008, 2013 

90-320 

2 

U Gem 

32 

1999, 2008 

16-320 

2,15 

CN Ori 

17 

1986 

137 

16 

IP Peg 

28 

1988, 2009 

50-150 

2,17 

CTCV J1300-3052 

24 

2010 

45 

18 

QZ Ser 

- 

2002 

200 

19 

Z Cha 

- 

1984 

150 

20 

HT Cas 

10 

2008 

60-190 

2 

OY Car 

28 

2010 

46 

21 

V2051 Oph 

31 

1998, 2009 

55-430 

2,23 

SDSS 103533.02+055158.3 

51 

2006 

70 

23 

WZ Sge 

384 

1996 

25 

24 

SDSS J143317.78+101123.3 

38 

2008 

32 

25 
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Table 2— Continued 


Object # Spectra 

Year 

Resolution (km/s) 

ref 

SDSS 1507+52 

2006 

160 

26 

El Psc = J2329+0628 

2001 

180 

27 


Note.- Information on SDSS 2044-04, RX 1951.7+3716, UY Pup, AT 
Ara, GY Hya, V392 Hya, RY Ser, SDSS 0813+45, LY UMa, TT Crt, CN 
Ori, QZ Ser, Z Cha, SDSS 1507+52 and El Psc has been extracted from 
literature. In the case of CN Ori, the FWHM and EW values are obtained 
from 17 individual measurements listed in table 1 of Barrera & Vogt (1989). 
The AH Her spectra are phased binned averages. 


Thorstensen (20051; (4) Watson et al. (2007); (5) Thorstensen et al.| ( |2004 |); 
(6) |Echevarria et al. (2007); (7) Bruch et al. (1997); (8) 


References. - (1) Martin et al. (1995); (2) this paper; (3) Peters & 


Casares et al. 


et al. 

(2001); (12 

Horne et al. (il986); (13) Pyrzas et al. (2012); (14) 

Southworth et al. 

(2009); (15) Echevarria et al. 

(2007); (16) Barrera & 

Vogt 

(1989); (17) 

Harlaftis et al. (1994); (18) Savoury et al. (20121; (19) 

Thorstensen et al. 

(2002a); (20) Marsh et al. (1987); (21) Copperwheat et 


al.|fl2012D; (22) |Steeghs etaT] (|2001D; (23) |Southworth et al.| fl2006P; (24) 


Skidmore et al. (2000), (25) Tulloch et al. 
^2008); (27) Thorstensen et al. (2002b) 


(2009); (26) Patterson et al. 
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Table 3 

X-RAY TRANSIENTS 


Object 

P(d) 

K 2 (km/s) 

FWHM (km/s) 

EW (A) 

Ml (Mg) 

i (deg) 

reft 


Black Holes 

V404 Cyg 

6.47129 

208.4 ± 0.6 

1029 ± 94 

19.0 ± 5.2 

Q f)+°- 2 
J ’ U -0.6 

67±? 

1,2 

BW Cir 

2.54451 

279.0 ± 4.7 

1062 ± 91 

53.4 ± 8.2 



3 

XTE J1550-564 

1.5420333 

363.1 ± 6.0 

1506 ± 151 

28.4 ± 5.8 

11.5±3.9 

75±4 

4 

N. Oph 77 

0.5228 

441.0 ± 6.0 

1885 ± 189 

68.5 ± 15.1 

6.2±1.2 

70±10 

5 

N. Mus 91 

0.4326058 

406.0 ± 7.0 

1735 ± 92 

48.7 ± 6.5 

7.0±0.6 

54±2 

6,7 

GS 2000+25 

0.3440915 

519.9 ± 5.1 

2207 ± 341 

22.4 ± 7.2 

5.5—8.8 

58-74 

8,9 

A0620-00 

0.32301405 

437.1 ± 2.0 

1966 ± 110 

36.4 ± 10.8 

6.6±0.3 

51±1 

10,11 

XTE J1650-500 

0.3205 

435.0 ± 30.0 

1866 ± 348 

8.7 ± 2.3 



12 

N Vel 93 

0.285206 

475.4 ± 5.9 

2082 ± 208 

33.7 ± 7.4 



13 

XTE J1859+226 

0.274 

541.0 ± 70.0 

2361 ± 109 

129.1± 17.8 



14 

GRO J0422+320 

0.2121600 

380.6 ± 6.5 

1716 ± 202 

165.3± 18.2 



15 

XTE J1118+480 

0.1699337 

708.8 ± 1.4 

2850 ± 187 

86.7 ± 18.9 

CN 

00 

1 

05 

CO 

68-79 

16,17 

Neutron Stars 

Cen X-4 

0.6290522 

144.6 ± 0.3 

678 ± 48 

36 ± 9 

1 04+0-37 
1 ' y ^-0.85 

32t® 

18,19 

XTE J2123-058 

0.24821 

298.5 ± 6.9 

1437 ±272 

18.5 ±5.0 

1.5±0.3 

73±4 

20,21 


t This column provides references for the adopted values of P, Jv 2 , Mi and i. 


References. — (1) Casares (1996); (2) Khargharia et al. (2010); (3) Casares et al. (2009a); (4) Orosz et al. (|2011 ); 
(5) |Harlaftis et ah] (1997); (6) 


et al. ( 

2004); (10) 

Gonzalez Hernandez & Casares (2010); (11) Cantrell et al. 

(2010); (12) Orosz et al.| (2004); (13) 

Filippenko et al. 1 

1999); (14) Corral-Santana et al. (2011); (15)|Webb et al. 

(2000); (16) Gonzalez Hernandez et 


al. (2012); (17) Khargharia et al. (2013); (18) Casares et al.| (2007); (19) |Shahbaz et al.| (2014); (20) Tomsick et al. 


(2001); (21) Zurita et al. (2000). 
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Table 4 

Cataclysmic Variables 


Object 

P(d) 

K -2 (km/s) 

FWHM (km/s) 

EW (A) 

ref t 

GK Per 

1.9968 

120.5+ 0.7 

722 ± 113 

18.1+ 4.0 

1 

SDSS 2044-04 

1.68 

90.0+ 8.0 

420 ± 29 

21.0 ± 2.9 

2 

BV Gen 

0.611179 

137.3+ 0.3 

820 ± 58 

13.0 ± 1.4 

3-4 

RX 1951.7+3716 

0.492 

81.0 ± 7.0 

566 ± 40 

29.0 ± 4.1 

2 

UY Pup 

0.479269 

102.0+ 4.0 

615 ± 43 

17.0 ± 2.4 

5 

EYCyg 

0.4593249 

54.0+ 2.0 

341 ± 63 

8.0 ± 3.0 

5 

DX And 

0.4405019 

105.8+ 3.8 

761 ± 96 

5.4 ± 0.9 

7-8 

SY Cnc 

0.3823753 

88.0+ 2.9 

560 ± 58 

10.9 ± 4.2 

9 

AT Ara 

0.37551 

99.5+ 3.2 

618 ± 43 

30.0 ± 4.2 

10 

RU Peg 

0.3746 

121.0+ 2.0 

881 ± 21 

21.4 ± 01.4 

11 

GY Hya 

0.3472309 

176.0+8.0 

1267+ 89 

9.0 ± 1.3 

2 

CH UMa 

0.3431843 

76.0+ 3.0 

518 ± 22 

46.3 ± 2.2 

5 

V392 Hya 

0.324952 

144.0+ 9.0 

850+ 60 

14.0 ± 2.0 

2 

RY Ser 

0.3009 

87.0+ 6.0 

519 ± 36 

22.0 ± 3.1 

5 

HS 0218+3229 

0.297229661 

162.4+ 1.4 

1023 ± 70 

27.9 ± 4.0 

12 

EM Cyg 

0.290909 

202.0+ 3.0 

1242 ± 86 

18.0 ± 0.9 

13-14 

Z Cam 

0.289840 

193.0+ 17.0 

1107 ± 64 

43.1 ± 1.3 

15 

SDSS 0813+45 

0.2890 

54.0+ 7.0 

461 ± 32 

30.0 ± 4.2 

5 

V426 Oph 

0.285314 

179.0+ 2.0 

1224 ± 13 

12.5 ± 0.3 

16-17 

ss Cyg 

0.27512973 

165.0+ 1.0 

986 ± 14 

39.5 ± 2.0 

17-18 

LY UMa 

0.271278 

141.0+ 3.0 

1002 ± 70 

28.0 ± 5.6 

19 

BF Eri 

0.270881 

182.5+ 0.9 

1069 ± 34 

27.7 ± 0.5 

20 

TT Crt 

0.2683522 

212.0+ 5.0 

1310 ± 92 

17.0 ± 2.4 

5 

AH Her 

0.258116 

175.0+ 2.0 

1037 ± 52 

16.2 ± 1.9 

18,21 

EX Dra 

0.20993698 

210.0+ 14.0 

1292 ± 63 

37.0 ± 4.0 

22-23 

HS 2325+8205 

0.194334535 

237.0+ 28.0 

1377 ± 93 

48.8 ± 7.2 

24 

DQ Her 

0.193620919 

227.0+ 10.0 

1168 ± 101 

28.8 ± 1.5 

25-26 

SDSS J100658+233724 

0.18591324 

258.0+ 12.0 

1446 ± 144 

91.7 ± 10.2 

27 

SS Aur 

0.1828 

167.0+ 15.0 

911 ± 30 

95.8 ± 5.2 

28-29 

U Gem 

0.17690619 

298.4 ± 9.0 

1602 ± 53 

87.7 ± 15.6 

30-31 

CN Ori 

0.16319006 

231.0+ 7.0 

1725 ±395 

23.8 ± 8.6 

31-32 

IP Peg 

0.1582061029 

331.3+ 5.8 

1899 ± 208 

99.8 ± 28.0 

33-34 

CTCV J1300-3052 

0.088940717 

378.0+ 6.0 

2064 ± 66 

127.2 ± 14.7 

35-36 

QZ Ser 

0.0831612 

207.0+ 5.0 

1218 ± 85 

36.0 ± 5.0 

37 

Z Cha 

0.0744992335 

398.0+ 9.0 

2000 ± 140 

79.3 ± 11.1 

38-39 

HT Cas 

0.0736472029 

389.0+ 9.0 

1912 ± 155 

154.4 ± 18.9 

40-41 

OY Car 

0.0631209343 

470.0+ 2.7 

2125 ± 147 

73.7 ± 9.0 

42-43 

V2051 Oph 

0.0624278634 

436.0+ 11.0 

1981 ± 107 

159.1 ± 14.2 

44-45 

SDSS 103533.02+055158.3 

0.0570067 

520.0+3.0 

1996 ± 125 

51.5 ± 8.3 

46 

WZ Sge 

0.056688 

510.0+ 15.0 

1986 ± 129 

150.1 ± 22.8 

47-48 

SDSS J143317.78+101123.3 

0.054240679 

511.0+2.0 

2460 ±148 

141.9+14.1 

46 
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Table 4— Continued 


Object 

P(d) 

K 2 (km/s) 

FWHM (km/s) 

EW (A) 

ref ^ 

SDSS 1507+52 

0.04625834 

529.0+17.0 

2325 ±163 

87.0+17.4 

49 

El Psc = J2329+0628 

0.044566 

348.0+ 4.0 

1865 ± 131 

25.0 ± 3.5 

50 


1 This column provides references for the adopted values of P and I\ 2 . 


References. — (1) Morales-Rueda et al. (2002); (2) Peters & Thorstensen (2005); (3) 
Gilliland (1982 ); (4)|Watson et al. (2007); (5) Thorstensen et al. ( 2004|): (6)~ Echevarria et al. 


(1 2007 I); (7) ]Bruchet al.|(|1997D ; (8)|Drewet al.|(|1993| ; (9)|Casares et al.|([2009b ) 


(2003); (11) Stover 


et al. (2000D; (15) 


,_ ( 10 ) 

(|1981D; (12)|Rodriguez-Gil et al.|fl2009j); ( 13) |Robinson| ([1974D ; (14) 
Kraft et al.|(1969|); (16) Hessman (1988 ); (17) |North et al. 


Hessman et al. (1984); (19) Tappert et al. (2001); (20) Neustroev & Zharikov 


Horne et al. (1986); (22) 


(2012); (25) Wood et al. 


Baptista et al. (2000); (23) Billington et al. (1996 ); (24) 
(2005); (26) Horne et al. ( 1993| ); (27)~ Southworth et al. 


Bruch 


North 


(2002); (18) 


(2008); (21) 


Pyrzas et al. 


(2009); (28) 


Shatter & Harkness (|1986); (29)|Davey &; Smith (19921; (30) Marsh et al.| (1990); (31) Friend 


et al. ( |1990| ); (32j Barrera & Vogt (1989 ); (33) |Copperwheat et al. ( 2010|); (34) Beekman et al. 


(2000); (35) Savoury et al.| (2011 ); (36) |Savoury et al. ~ p012 ); (37) Thorstensen et al.J j 2002a|) ; 


(38) |Wood et al.|(|1986D; (39) |Wood k Horne l (|199(f; (40) |Borges et al | Q2008| ; (41) |Horne 
et al. (1991); (42) Greenhill et al. ( |2006 ); (43) Copperwheat et akf (2012); (44) Baptista et al. 


(2003); (45) Baptista et al. (1998 ); (46) |Littlefair et ah ( |2008 ); (47) Patterson et al. |( 19981; 
(48) |Steeghs et al.' (2007); (49) Patterson et al. (2008); (50 jThorstensen et al. (2002b). 
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